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Quenched SU(3) hadron spectroscopy using improved fermionic and 
gauge actions 

Sara Collins, Robert G. Edwards*, Urs M. Heller, and John Sloan 
SCRI, Florida State University, Tallahassee, FL 32306-4052, USA 

We present results of quenched SU(3) hadron spectroscopy using O(a) improved Wilson fermions. The config- 
urations were generated using an 0(a 2 ) improved 6-link SU(3) pure gauge action at /3's corresponding to lattice 
spacings of 0.43, 0.25, 0.20, 0.18, and 0.15 fm. We find evidence that fermionic scaling violations are consistent 
with 0(a 2 ) errors. 



1. INTRODUCTION 

The Symanzik action improvement program 
has been proposed |l],|2| as a way to reduce scal- 
ing violations in the approach to the continuum 
limit from a lattice action. In another contri- 
bution, we report on our preliminary investiga- 
tions into the nature of scaling violations inher- 
ent in improved actions. In another proceed- 
ings in this volume pj, we make detailed com- 
parisons of Wilson and tadpole-improved Clover 
fermion actions on configurations with dynamical 
fermion content. In this work, we perform SU(3) 
quenched spectroscopy using a one-loop tadpole- 
improved Symanzik pure gauge action and the 
tree-level tadpole-improved Clover fermion ac- 
tion. We measure the hadron spectrum and the 
string tension at lattice spacings 0.15, 0.18, 0.20, 
0.25, and 0.43 fm. The goal of this work is to mea- 
sure the lattice spacing dependence of the hadron 
spectrum. We find significant scaling violations in 
arciK* I V a 2 <J consistent with Oia 2 ). 

2. ACTIONS 

Recently, it was shown [|| that the tadpole im- 
proved Symanzik inspired pure gauge action Q is 
successful in restoring Lorentz invariance in the 
static quark potential even at a lattice spacing 
of 0.43fm. Inspired by this result, we used the 
action 
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= -(l + 0.4805a s ) / 20ul 
= -0.03325cv s / ul . 
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The sums are over plaquettes (pi) ,1x2 rectangles 
(rt), and the six length path (/z, v, p, — p, —v, —p) 
with p,v,p all different (pg). We use the mea- 
sured expectation value of the plaquette to deter- 
mine ito- The coupling constant a s is calculated 
using the tree level action, 
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3.06839 

Classically, this action has 0(a A ) errors but 
quantum effects induce 0{a 2 a 2 ) errors. However, 
it was shown [|| that the action (|l|) is insensitive 
to nonperturbative tuning of the coefficients sug- 
gesting that the 0(a 2 a 2 ) errors are small. 

To reduce fermionic scaling errors, we use the 
tadpole-improved version of the tree-level Clover 
fermionic action proposed in 

qSW 



Op 
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where is the usual Wilson fermionic action 
and F p i is a lattice definition of the field strength 
tensor. For the tree level action, c(g 2 ) = 1. 
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Figure 1. Scaling plot of am,K* versus the string 
tension a 2 a. The horizontal line is the physical 
value using a 2 a = 440MeV. The diagonal line is 
a least squares fit. 
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Figure 2. Scaling plot of Jd versus the string 
tension a 2 a. The physical value is 0.499. 



With this tree-level action, one expects quan- 
tum errors of 0{aa). However, Naik |^| computed 
the one loop and the dominant two loop contri- 
butions to c((?q) and finds that after tadpole im- 
provement the coefficient of g\ is 0.016. There- 
fore, we expect to see only the dominant 0(a 2 ) 
fermionic errors. 

3. OBSERVABLES 

To set the scale, a, we used the string ten- 
sion. We computed finite T approximations to 
the static quark potential using time-like Wil- 
son loops W(R,T) which were constructed us- 
ing APE'-smeared spatial links ||. On and 
off-axis spatial paths were used with distances 
R = n, V2n, V3n and y/bn, where n is a posi- 
tive integer. The "effective" potentials were fitted 
using a correlated \ 2 procedure with the spatial 
covariance estimated using bootstrap. The string 
tension was estimated using the form 



V{R) = V + <jR 
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Figure 3. Scaling plot of the nucleon versus the 
rho as a function of the string tension a 2 a. The 
horizontal line is the physical value of 1.22. 
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The last term takes account of the lattice arti- 
facts at short distance. Here, Gl(R) denotes the 
lattice Coulomb potential for the Wilson gluonic 
propagator. We did not have available the lattice 
Coulomb potential for the tree-level Symanzik ac- 
tion; however, we found Gl(R) enabled fits to 
smaller R. 

For hadron measurements, we used correlated 
multi-state fits to multiple correlation functions 
as discussed in [Q . We computed quark propaga- 
tors in Coulomb gauge at four n values for each (3. 
We used two gaussian source smearing functions 
with smeared and local sinks. The widths were 
chosen to have a positive and negative overlap, 
respectively, with the first excited state. This im- 
proved the signal for the excited state in a multi- 
state ("vector") fit. We also measured mixed va- 
lence content mesons and baryons. 

As discussed in Q , we tried to avoid human in- 
tervention as much as possible in the fitting proce- 
dure. For the multistate fits, we chose the plateau 
with max NdofQ J\ i 6 Mi /Mi where the product is 
over the relative errors of the ground and first ex- 
cited state masses, and Q is the confidence level 
for a Ndof degree of freedom fit. 

Once a fitting range was chosen, a bootstrap 
ensemble of masses was generated and used for 
correlated \ 2 chiral fits. We fit the particle mass 
according to the ansatze 
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This ansatz was motivated by one loop chiral per- 
turbation corrections to baryon masses. We note 
that recently one loop Mp s corrections to the vec- 
tor mass were derived in [Q. 

An important method used to increase statis- 
tics is via 'Z(3)' fermion sources ||. Within the 
source time-slice source for quark propagators, 
smearing functions are placed at regular inter- 
vals along each spatial dimension. These smear- 
ing functions are weighted with random elements 
chosen according to a Z(3) distribution. Since 
the elements are not correlated, each origin in a 
mesonic or baryonic construction contribute in- 
dependently to the correlation functions. The 
distance between origins should be chosen larger 
than the lattice size of a physical state or the 



smearing width used. 

4. SIMULATIONS 

We generated quenched configurations using 
the action (Q) at = 7.90, 7.75, 7.60, 7.40, and 
6.80. Unfortunately, the non-local nature of the 
action precluded an efficient implementation of an 
over-relaxed / heat-bath algorithm for the CM-2 
at SCRI - we would need to divide the lattice into 
32 sublattices. We used a Hybrid Monte Carlo al- 
gorithm with a second order integration scheme, 
step size 5 = 0.1, and a global acceptance test. 
We typically found > 0.75% acceptance rates. 

We determined meson and baryon masses us- 
ing the action (g) . The pseudoscalar masses were 
tuned to the same physical mass (in units of 
V a 2 er) between the various (3. The vector to pseu- 
doscalar mass ratios varied from about 2.0 down 
to 1.2. At the coarsest lattice, we found a large 
number of exceptional configurations on an 8 3 x 16 
lattice. To avoid problems, we used a 16 3 x 32 lat- 
tice throughout except at (3 = 7.40. We expect 
finite volume effects to be small since the box size 
is 2.40fm at (3 = 7.90. 

The results of the chiral extrapolations of our 
best fits are listed in Table [l] along with the string 
tension. Note the lattices are too coarse to reli- 
ably use Sommer's force method since the lattice 
spacing is close to the empirical 0.5fm scale. The 
lattice spacings we quote used V a 2 a = 440MeV. 
For the vector mesons, we used a quadratic fit 
ansatz of (§) for (3 = 7.90, 7.75, and 7.60 using 
all 10 mixed valence states. At (3 = 7.40, we 
dropped the two largest masses and used a for a 
quadratic fit. 

At (3 — 6.80, we found some interesting re- 
sults. First, the assumption of using mixed va- 
lence states as distinct measurements appeared to 
break down - the meson masses fluctuated signif- 
icantly enough to lower the quality of the fit. For 
this reason, we only used the four diagonal masses 
(same valence content) for fits. Also, the smallest 
pseudoscalar measured had amps = 0.9939(13). 
At these large masses, we expect significant devi- 
ations from the continuum dispersion relations. 
Quite significantly, a quadratic fit ansatz was 
clearly violated, and we used a cubic fit ansatz 
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Table 1 

String tension and best fits of masses extrapolated to physical ratios. 



p 


L 


■^Vi meas 


MO 


a 1 a 


am p 


amK* 


Jd 


am,N 


7.90 


16 


100 


0.8848 


0.1272( 16) 


0.6206( 58) 


0.6937(44) 


0.3724(58) 


0.881( 1) 


7.75 


16 


100 


0.8800 


0.1658( 50) 


0.7123( 46) 


0.7950(35) 


0.3679(40) 


0.869(25) 


7.60 


16 


32 


0.8736 


0.2184( 13) 


0.7666( 85) 


0.8565(68) 


0.3707(71) 


0.888(46) 


7.40 


8 


68 


0.8629 


0.3470( 29) 


0.9482(115) 


1.0577(89) 


0.3658(75) 


1.154(59) 


6.80 


16 


72 


0.8261 


0.8961(177) 


1.2706( 71) 


1.4120(41) 


0.3543(30) 


1.778(70) 



of (||) with Ndof = 1. This indicates that scal- 
ing violations inherent with large am q are roughly 
approximated by higher order terms in the chiral 
ansatz of (||). 

In our chiral fits for the nucleon, we only used 
the four diagonal masses. For all propagator fits, 
we used a two exponential fit with two smearings. 
At j3 = 7.90 - 7.40, we found that the quadratic 
ansatz broke down at the largest mass. However, 
there was not sufficient signal for a cubic term. 
At j3 = 6.80, we used a cubic ansatz for all four 
masses. At j3 = 7.90, the extrapolated nucleon 
value is way off scaling. We think this comes from 
a lack of statistics due to autocorrelations, and 
the error bar is underestimated. 

5. SCALING VIOLATIONS 

Our main result for the determination of scal- 
ing violations is in Figure [l]. We interpolate (or 
extrapolate) the vector meson mass to the K* /K 
ratio and find that there are significant scaling 
violations in the arriK* /V a 2 a ratio. Fitting to a 
line, we find Q — 0.1 and the intercept of 1.99(2) 
which is close to the physical value of 2.03 us- 
ing \J a 2 a — 440MeV. This implies that the as- 
sumption of 0(a 2 ) errors is significant. If we 
fit to V a 2 a assuming 0(a) errors, we find that 
Q = 0.17 but the intercept is 2.22(3). This value 
would imply \J a 2 a = 347MeV which is probably 
too low. Indeed, constraining the 0(a) fit to pass 
through the value 2.0275 has Q = 10~ 9 . We thus 
feel reasonably confident that the fermionic scal- 
ing errors are (at least) consistent with 0(a 2 ). 

Following M, we plot the ratio 



versus a 2 a in Figure ||. We see that Jd is reason- 
ably scaling out to j3 = 7.40 indicating that scal- 
ing violations in the individual masses are nearly 
canceling in the ratio. The value of Jd is consis- 
tent with other quenched calculations. 

In Figure [| we plot the ratio of arriN /am p . We 
see reasonable scaling out to f3 — 7.40, but the 
value at (3 = 7.90 is way off as described above. 

The computations done on the CM-2 at 
S.C.R.I. This research was supported by DOE 
contracts DE-FG05-85ER250000 and DE-FG05- 
92ER40742. 
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